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ABSTRACT 


Data obtained with the Mariner 9 infrared spectroscopy experiment during 
the Martian Dust Storm of 1971-72 are examined for evidence of topographic 
influence on the atmospheric thermal structure. Temperature perturbations 
which are well correlated with the planetary scale topography are found super- 
posed on the large amplitude diurnal thermal tide previously reported. Appli- 
cations of tidal theory to the data indicate that the observed perturbations result 
from the Idnematic interaction of the westward traveling diurnal wave with the 
large scale components of topography. The dominant mode is excited by the 
wave-number two topography component and is a vertically evanescent eastward 
traveling wave with an equivalent depth comparable to the atmospheric scale 
height. The principle dynamic effect of this mode is the enhancement of the 
amplitude of the near-surface diurnal wind to over 40 m s“^ in limited areas 
near 30° south latitude. Thus, it appears likely that dust was injected 
into the atmosphere in these regions during the storm. Other waves excited 
include vertically propagating modes which produce little effect in the lower 
atmosphere, but may represent a source of energy for the upper atmosphere 


during dust storm conditions. 


Introduction 


The Mariner 9 spacecraft was injected into orbit about Mars on 14 Novem- 
ber 1971 during a major planet-wide dust storm. Extensive data were acquired 
during the storm and its subsequent decay. Analyses of these data have provided 
information on the physical and dynamical state of the atmosphere during the 
storm. The dust cloud which extended to a height of ~50 km (Leovy, et al. 1972) 
was found to exert a profound influence on the thermal structure of the atmos- 
phere (Hanel, et al., 1972; Kliore, et al., 1972). Atmospheric temperatures ob- 
tained with the Mariner 9 infrared spectroscopy experiment (IRIS) indicated 
mean temperatures substantially higher than those expected for a clear COj 
atmosphere. A diurnal temperature wave was fmmd with a peak-to-pealc ampli- 
tude approaching 30K at some latitudes extending from the surface to at least 
40km (Hanel, et al., 1972; Conrath, et al., 1973). 

Temperature data obtained during the dust storm have been used to study 
the thermally driven atmospheric tidal system (Leovy, et al., 1973; Pirraglia 
and Ct math, 1974). In these studies, the surface of the planet was assumed to 
be a smooth sphere. However, it is known from a variety of measurements that 
the surface of Mars possesses substantial topographic relief on both a local and 
planetary scale (see for example Downs, et al., 1975; Kliore, et al., 1972; Hord 
et al., 1974; Conrath, et al,, 1973). This relief can locally reach or exceed an 
atmospheric scale height. Lindzen (1970) has pointed out that the atmospheric 
tidal waves can interact with the large scale topography to produce perturbations 
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on the Martian tidal regime, and Zurek (1976) has considered this problem in 
some detail. It is the purpose of the present study to examine the Mariner 9 


IRIS temperature data for evidence of the existence of effects related to the 
large scale topographic relief. 

First the IRIS temperature data will be analyzed and correlated with topo- 
graphy. Next a first order model of the modulation of the tidal waves by large 
scale topography will be developed. Finally the theoretical results will be com- 
pared with the observational data, and conclusions will be drawn concerning the 
effectiveness of the large scale topography in modifying the tidal structure. 

Data Analysis 

Temperature profiles were reconstructed from the IRIS data using meas- 
urements in selected spectral inteiwals within the 667 cm~^ COj absorption 
band. The analytical methods employed are similar to those developed for use 
in the terrestrial meteorological satellite program (see for example Conrath 
and Revah, 1972). The vertical resolution of the retrieved profiles is about 
one-half of an atmospheric scale height and the RMS random error due to in- 
strumental noise is estimated to be 2-3K. Useful information is obtained from 
the planetary surface up to about the 0.1 mb level. 

Temperature at the 2 mb level of the atmosphere is plotted versus longitude 
for three different latitude intervals in Figures 1,2, and 3. Each data point is 
talten from an individual temperature sounding. The data were obtained during 
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the first 47 days of the Mariner 9 mission and are shown for both morning and 
evening conditions. At the bottom of each figure, the mean surface elevation 
is given for the latitude interval 0° to -45° for purposes of comparison. These 
data were also obtained from the IRIS experiment using measurements in the 
wing of the 667 cm'^ CO^ band (J. Pearl private communication, 1975). In 
both the early morning and late afternoon, there is an apparent correlation with 
topography. The atmosphere is warmer over depressions and cooler over higher 
elevations in the morning. In the late afternoon, the situation is reversed with 
warmer temperatures occuring over higher elevations and cooler temperatures 
over lower regions. The mid-day data coverage is much less complete and 
a detailed analysis cannot be made during the local time intervals from mid- 
morning to early afternoon. 

One possibility for the observed temperature behavior is the modulation 
of the thermal tide by tlie large scale topography. To persue this possibility 
further, it is necessary to first consider the theory of topographic interaction 
with the tidal wave. 

Theory 

Chapman and Lindzen (1970) have given a qualitative discussion of the 
general problem of the interaction of tides with the planetary surface, demon- 
strating that waves may occur with phase speeds different from that of the solar 
excitation, and Zurek (1976) has developed a computational metliod for treating 
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the complete tidal problem. For purposes of Interpreting the data considered 
here, however, a simplified theoretical model wll be developed. The approach 
talcen is as follows. The amplitude of the diurnal temperature wave can be re- 
garded as consisting of a zonal mean plus a longitudinally varying component. 

The zonal mean will be used, along with the observed planetary- sc ale topography, 
to calculate the longitude dependent component. A comparison can then be made 
with the observed longitudinal variations. 

The theory of the thermally driven atmospheric tide on a smooth spherical 
planet has been reviewed in detail by Chapman and Lindzen; for a discussion of 
the basic assumptions involved and the detailed mathematical development, the 
reader is referred to their work. The basic independent variable normally 
employed in atmospheric tidal theory is 


1 ^ 
^Po 


( 1 ) 


where p is atmospheric pressure, Po is atmospheric pressure for the unper- 
turbed (static) state, and 7 is the ratio of specific heats. If it is assumed that 
J, the thermotidal heating per unit mass per unit time, can be written 


J(^^. <P, z. 


s 


t) = J ' ° ( 8 , Z) e 


> (Sg^ + CTt ) 


(21 


and that 


G(e, 0 , Z, t) = 0, Z) 


(3) 


4 


7 


then it can be sliovvn that G°^ must satisfy 
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The operator is defined as 
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The notation is as follows; 
d - colatitude 
0 = east longitude 
a = wave frequency 

Sq = longitudinal wavenumber of source of excitation 
Z = vertical height 
H = scale height 
a = radius of planet (=3383 km) 

n = siderial rotation rate of planet (=7.09003 X 10"® sec"*) 
g = gravitational acceleration (=370 cm sec"^) 

« = R/Cp (=0.23) 

R = gas constant (=1.89 x 10"® ergs gm“* deg"*) 

Cp = specific heat at constant pressure (=8.22 x lo"® ergs gm 
f = a/212 


(4) 
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For a diurnal wave with a period of one mean solar day, a - 7.07765 x 10 
sec"^ . The scale height H will be tal^en as constant in all subsequent calculations, 
and the values used for the thei’modynamic parameters are for a pure CO^ 
atmosphere. 

With the specification of appropriate boundary conditions, (4) can be 
solved for G from which the other tidal field variables can be obtained. In 
addition to the requirement that the solution be finite at the poles, boundary 
conditions must be specified at the planetary surface and for z - oo . The upper 
boundary condition will be discussed later. The lower boundary condition for 
a smooth spherical planet is usually taken as w (z = 0) = 0 where w is the vertical 
tidal velocity. In the presence of topography, this condition must be generalized 
to the requirement that the component of tidal velocity locally normal to the 
surface must vanish. If ^ {B, cj)) represents the surface elevation above a 
reference surface of constant geopotential, this condition can be written approxi- 
mately as 


w = V • a t z = 0 ^0) 

where v is the horizontal velocity field and V is the horizontal gradient operator. 

Consider now a component of topography with longitudinal wave number m 
written in the form 


<p) - ^JO) cos(m0 + 


( 7 ) 
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Let the following definitions be intl.Tduced; 
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b_ (6<) = bV^^) 

where the asterisks denote complex conjugation, and e = max /H. Introduc- 
tion of (7) into (6) gives the appropriate lower boundary condition 






( 8 ) 


where u is the southward component of tidal velocity and v is the eastward 
component. 

If it is assumed the a^^i and b^^ are of order unity and »^ < < 1, then a sep- 
arable problem can be obtained by expanding G"^, u'^ , and in powers of , 
i.e., 
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ir = 11^ + eu^ + + ... 


( 10 ) 
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Substituting (8) into (4) and equating like powers of e yields an infinite set of 
differential equations. The zero order and first order equations are 
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and 
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where x = z/M. In a similar fashion, substitution of (10) and (11) into (8) gives 
the correspc nding lower boundary conditions 
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Considering first the zero order tidal equation (12), separation of variables is 
accomplished by assuming 
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wdiere the functions satisfy the eigenvalue equation 
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and the height dependent factors are obtained from the vertical structure 
equation 
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Equation (18) is Laplace's Tidal Equation, and the eigenfunctions 0]^' ^ are 
Icnown as Hough functions. They can be represented as infinite series of as- 
sociated Legendre functions, and their properties have been extensively dis- 
cussed in the literature. The eigenvalues are here written in the conventional 
form -4a^H^/gh^j'^° where h^’^° is the so-called equivalent depth. 

The tidal fields which are of interest for the present discussion can be 
written 
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where the superscripts a ,s have been surpressed. and ^ p are the tidal 
temperature and pressure fluctuations. The functions and V, are related 
to by 


U = i +T^cot tJ H ( 0 ), 

- cos^ fi ^ / " 


M ^ I cos d ^0 1 , , ^ 

f (26) 

The height dependent factors in (20) -(24) can be expressed in terms of solutions 
of the vertical structure equation (19): 
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Using (31), the zero order lower boundary condition (14) can be written 



dx 


+ 




at X = 0. 


(32) 


In the usual approach to the tidal problem, the J^' s are specified, and (19) is 
solved subject to (32) ,ond an upper boundary condition. However, in the present 
problem, the temperature field can be directly specified from the observations. 
It is desired to talte the zero order tidal temperature field as given, 
and calculate the perturbation fields due to interaction of the zero order 
fields with the topography. The zero order quantities which will be 
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reqvlired in obtaining the first order solution are (0) and Vj^(O), and from 
(29) and (32) 


4an2 


A relationship between y^'^°(0) and<^T^'^°^ , the density weighted vertical 

Sn 

average of the zero order tidal temperature component ST„ i is derived in 
the Appendix. Using this result, (33) can be wi’itten 
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Solution of the first order equation (13) subject to the lower boundary con- 
dition (15) will now be considered. Let 


G‘[(6'. 0, x) = x) 


and substitute into (13). This yields a set of homogeneous equations of the 


form 
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where P_ is the differential operator in square brackets in the left hand side 
of (18) with So replaced by s. If G*^’ ^ is expanded in terms of where 
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there results a set of homogeneous vertical structure equations 


dx2 4 


The first order vertical velocity can be expanded in the form 




Using this expression along with the appropriate expansions for Vg and Uq 
and substituting into the first order lower boiindary condition (15), it can be 
seen that the only non-zero terms in the sum over s in (39) correspond to 
s = Sg ± m. MaJdng use of the orthogonality properties of the Hough functions, 
the boundary condition becomes 


■•77 j 
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In obtaining this expression, use has been made of (30). The right hand side 
of (40) can be expressed directly in terms of the vertical mean zero order 


tidal temperature field using (34) ; 
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Explicit solutions of the vertical structure equation (38) with the inhomo 
geneous lower boundary condition (41) must now be considered. With the 
definition 



the general solution can be written 



where and are arbitrary constants. For the case when k H/h - 

1/4 > 0 (vertically propagating solutions), it will be assumed that no downward 
propagation of energy occurs as x - “ with the result that ;= 0. When 
K H/h - 1/4 < 0 (vertically evanescent solutions), the requirement that the 
kinetic energy density remain bounded as x “ again leads to B^’^ =0, 
can be evaluated by substitution into (41), using (31) to express the vertical 
velocity in terms of The result Is 
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From (28), the height dependent coefficients for the first order tidal temperature 
field can be written 
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Finally, the total first order tidal temperature field produced by the interaction 
of the zero order tidal wave of wavenumber Sq with topography of wavenumber 
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The first order tidal surface pressure and velocity fields can be written 
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From (47) -(50) it can be seen that the first order effect of the presence of large 
scale topography is to excite waves with longitudinal wavenumbers equal to 
both the sum and the difference of the zero order wavenumber and the topo- 
graphy wavenumber. The phase speeds of these waves will in general be differ- 
ent from the speed of the solar excitation. Note that for m > s^, eastward 
propagating as well as westward propagating waves can occur. 

The interaction of the tidal regime with large scale planetary topography 
has been studied here by considering the Idnematic effects of the lower boundary. 
The calculations have been carried through first order in the ratio of the surface 
elevation to the atmospheric scale height. For Mars, this ratio, while less than 
unity on a planetary scale, is not particularly small. However, the model should 
serve as a useful tool for achieving an xmderstanding of the observational data. 


Comparison of Theoretical and Obsei’vational Results 

The theoretical derivations of the preceeding section will now be applied 
to the Mariner 9 observations. The zero order diurnal temperature amplitude 
and phase is approximated by a model derived from the Mariner 9 IRIS data 
by Pirraglia and Conrath (1974), Its representation in terms of Hough functions 
is shown in Figure 4; the individual Hough amplitudes and phases are listed in 
Table 1, The phase has been chosen such that tt/ 2 corresponds to a diurnal 
wa.ve maximum at 18 hours local time. The Hough function notation corresponds 
to that employed by Chapman and Lindzen (1970), Negative subscripts denote 
modes corresponding to negative equivalent depths. 

The large scale topography is approximated by the truncated expansion 


^/H = 



(C^ cos ni0 + sin m0) P^(^) 


(54) 


where P^(^) is the (unnormalized) associated Legendre function as defined by 
Jalmke and Emde (1945). Longitudinally symmetric (m = 0) components of 
topography can modify the migrating diurnal tidal wave, but cannot excite new 
tidal modes through first order interaction. The effects of the m = 0 compon- 
ents cannot be removed from the zonally averaged data on which the model of 
Table 1 and Figure 4 is based. Thus the "zero order" temperature field used 
in the calculations already includes the contributions from the zonally symmetric 
topography, so these terms have been omitted from equation (54). The coeffi- 
cients used in the expansion were tal<en from preliminary harmonic fits to the 
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geoid and topography (E, J. Christensen, private oomm unication, 1974) and 
are given in Table 2. 

The first order perturbation tidal fields were calculated for each of the 
six components of topography contained in (54). In each case the eight lowest 
order Hough functions were retained in the expansion in (45). The total first 
order temperature perturbation was calculated by summing the contributions 
from all topography components. The results for the 2mb level are shown as 
solid curves in Figures 1, 2, and 3. A surface pressure of 5mb and a scale 
height of 11km were assumed. As plotted in the figures, the first order per- 
turbations are combined ^vith the zero order diurnal wave so a direct com- 
parison with the observational data points can be made. The amplitude and 
shape of the calculated perturbations appear to be in reasonable agreement 
with the observations. 

While a total of forty-eight Hough modes are included in the calculations 
shown, only a few mal^e significant contributions. Those modes with amplitudes 
at 2mb in excess of 1°K are listed in Table 3. The 0 contribution which 
dominates the solution, is a wave nmnber one mode with a phase velocity equal 
in magnitude but opposite in direction to the solar excitation. This mode 
corresponds to a positive equivalent depth, but with h > 4 k h, and is therefore 
vertically evanescent. The phases of the tidal fields associated with this type 
of mode are independent of height while the amplitude grows exponentially. 
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Examination of temperature data up to the upper limit of the IRIS sounding 
range at ~0,1 mb tends to confirm this behavior. However, at the upper levels 
the calculated temperature perturbation amplitude tends to exceed that indi- 
cated by the data. This is possibly due to the omission of dissipative effects 
from the theoretical model v/hich would tend to limit the amplitude at high 
levels. The remaining three modes listed in Table 3 are also eastward travel- 
ing waves, and since their equivalent depths satisfy the condition 0 < h < 4 kH 
they are vertically propagating. The contributions of these waves to the total 
first order perturbation fields in the lower atmosphere are small compared to 
the 0 mode. However, they are of some interest because they are poten- 
tially capable of transporting energy from the lower to the upper atmosphere. 

The strong excitation of the 0 j-“* mode is related to the fact that the 
equivalent depth of this mode is of the same order of magnitude as the atmos- 
1 'lierlc scale height. Examination of (45) along with (43) indicates that in the 
absence of dissipation, resonance can occur when H = (1 -k) h. For k = 0.23, 
the resonant scale height would be 9.44 km corresponding to a mean basic state 
temperature of 185 °K. In the present model, the value H = 11 km (corresponding 
to a basic state temperature of 215®K) is sufficiently far from the actual reso- 
nance point that the neglect of dissipative effects probably does not introduce 
serious error in the calculation of the coefficient A in equation (45). For 
example an estimate assuming linear radiative damping indicates that A changes 
by less than 10 per cent for a damping time of 1 day. 
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Examination of the matrix elements , which provide measures of the 
relative efficiency of the various zero order perturbation modes for exciting 
1st order modes, indicates that the interaction of the * zero order mode 
with the component of topography is primarily responsible for exciting the 
mode. represents the lowest order component of topography with 
longitudinal wave number 2 and is symmetric with respect to the equator. 

' is the lowest order diurnal tidal mode with positive equivalent depth and 
is also symmetric with respect to the equator. The strong excitation of the 
0j'"i mode by the wavenumber 2 topography on Mars was previously predicted 
by Zurek (1976). 

Having identified the component of topography as the principle source 
of excitation, is nov/ possible to investigate the effects of this excitation on the 
atmospheric dynamics during the dust storm conditions. The amplitudes and 
phases of the near-surface tidal wind components resulting from the p\ topo- 
graphy term, as calculated from (49), (50), (52), and (53), are shown in Figures 
5 and 6. The zonal wind component of this eastward traveling 1st order wave 
is dominant at low latitudes, reaching a maximum amplitude of about 16m sec"* 
at tlie equator. Near-surface wind amplitudes and phases of the zero order 
diurnal tidal wave, based on the model of Pirraglia and Conrath (1974), are also 
shown in the figures. The combination of the westward traveling zero order 
wave with the eastward traveling first order wave results in a diurnally fluc- 
tuating wind field, the amplitude of which is dependent on longitude as well as 
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latitude. This dependence is displayed in Figures 7 and 8. The presence of 
the first order wave enhances the diurnal amplitude in some locations wliile 
diminishing it in others. Of particular interest is the Increase in zonal wind 
to amplitudes greater than 40m sec"' over relatively small areas near -30“ 
latitude. The tidal winds will be augmented by the presence of a zonally 
symmetric wind component (Pirraglia 1975) and possibly by topographic effects 
on a local scale. This suggests that the free stream threshold velocity for the 
injection of dust into the atmosphere may have been exceeded in these regions 
during some time period each day throughout the coiarse of the storm. 

Returning now to the three vertically propagating waves listed In Table 3, 
it is of interest to calculate the associated vertical energy flux. The time 
averaged flux associated with each mode can be written (Siebert, 1961) 


W) = 


T^b^PoCO) 
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(55) 


The combined energy flux of the three vertically propagating modes is shown 
as a function of latitude in Figure 9. For comparison, the vertical energy flux 
in the earth's atmosphere due to the diurnal tide is estimated by Lindzen (1967) 
to be ~ 7 ergs seo~'^ cm"2 at the equator. The flux shown in Figure 9 repre- 
sents the energy propagating upward from the surface which serves as a source 
of excitation for these modes. Tlie extent to which the energy actually propagates 
into the upper atmosphere is an open question at this point. To study this aspect 
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further would require a more detailed modeling of the higher levels. However, 
the calculations presented here indicated that topographically excited vertically 
propagating modes do represent a potential source of energy for the upper at- 
mosphere during dust storm conditions. 

Summaiy and Conclusions 

The Mariner 9 HUS data during the 1971 dust storm show well defined 
topography-correlated temperature perturbations superposed on the strong 
diurnal thermal tide. These perturbations appear to result from the topographic 
excitation of additional tidal modes which do not follow the sun. Using a first 
order perturbation calculation, it was found that the topography component of 
longitudinal wave number 2 provides a strong source of excitation. In particular, 
the Wj'"' mode was found to be dominant. This is an eastward traveling wave- 
number one mode witl^ an equivalent depth comparable to the atmospheric scale 
height. The good quantitative agreement between the theoretically calculated 
and observed temperature perturbations indicates that the principal tidal modes 
have been properly identified. The results presented here support the conclusion 
of Zurek (1976) that the mode should be strongly excited. 

Estimates of the near surface tidal wind fields indicate strong contributions 
from topographically excited modes at low latitudes. Near 30° south where 
maximum wnds are found, the amplitude of the dlurnally fluctuating zonal \vind 
is found to be increased to over 40 m sec"’ in limited areas. Although the 
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quantitative results obtained here must be regarded with due caution because 
of the crudeness of the model, they do Indicate tliat the presence of large scale 
topography may enhance ^vlnds sufficiently at some locations to contribute to 
the Injection of dust into the atmosphere during dust storm conditions. 

The calculations indicate that vertically propagating modes are also excited 
by topography. Wliile these modes produce little significant effect on the tem- 
perature and wind fields of the lower atmosphere, they do result in a strong 
energy flux upward from the surface. These modes thus represent a possible 
source of energy for the upper atmosphere during a storm. However, to esti- 
mate the extent to wiiich this energy actually propagates to the upper levels 
requires a more detailed modeling of the atmosphere than has been considered 
hero. 
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APPENDIX 


The vertical structure equation for an isothermal atmosphere (19) and the 
expression for the corresponding tidal temperature field (28) can be combined 
to eliminate J yielding 


- i y ^ 

dx2 4 7Kh 

where all subscripts and superscripts designating modes have been surpressed. 
This can be regarded as a vertical structure equation with T(x) as the driving 
term. An expression for , the density weighted average of ^ T, can be 
obtained by multiplying both sides of (A-1) by and integrating to obtain 

<»T> = f <■-«« i!!>: clx -1 C c ydx. ,a-2) 
Jo dx2 4 Jo 

Integrating the first term on the right hand side by parts twice and assuming 






» 0 


X-Oj 


gives 
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(A-3) 


Tlie zero order lower boundary condition (32) can be employed to eliminate 
(dy/dx)^^o resulting in 
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TABLE 1 

Temperature amplitudes and phases for the Hough function representation 
of the zero order depth averaged diurnal thermal tide. 


Hough 

Mode 

Equivalent 

Depth 

(KM) 

Amplitude 

(K) 

Phase 

01.1 

.49 

3.60 

- '/t/2 


.085 

.94 

77/2 

01.1 

1044.2 

4.40 

7t/2 

Ql.l 

-8.66 

11.44 

77/2 

01,1 

-5.02 

3.46 

- 77 /2 


TABLE 2 

Coefficients for spherical harmonic representation of the ratio of 
the surface elevation to the atmospheric scale height ^/H. 


m 

't 



1 

1 

.00326 

-.13786 

1 

2 

.04050 

-.00620 

1 

3 

-.00242 

-.04287 

2 

2 

-.06410 

.04943 

2 

3 

-.01072 

.01249 

3 

3 

.00545 

.00910 
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TABLE 3 

The Hough modes most highly excited by topography. The topography 
mode responsible for the excitation is given in the final column. 


Hough 

Mode 

Temperature 
Amplitude at 
2mb (IC) 

Phase 

(Radians) 

Equivalent 
Depth (KM) 

Vertical 
Wave Length 
(Scale Heights) 

Topography 

Mode 

e|.-i 

7.09 

-2.23 

12.3 

3.45* 

Pa 

yt.-t 

1.11 

-2.89 

0.369 

2.44 

Pa 


2.01 

3.04 

3.45 

8.98 

p3 

01.-2 

1.30 

-0.87 

0.295 

2.18 

P 3 


Vertically evanescent. The volue given is the vertical distance over which the amplitude 
changes by a factor of e. 














FIGURE CAPTIONS 


Figure 1 


Figure 2 
Figure 3 
Figure 4 

Figure 5 
Figure 6 


Comparison of atmospheric temperatures at the 2mb level with 
topographic relief. The temperature data points represent in- 
dividual measurements obtained within the latitude interval -20° 
to -30° for the local time intervals indicated. The solid lines 
represent the results of model calculations for -25° latitude and 
6.5 and 17.5 hours local time. The mean surface elevation ex- 
pressed in units of atmospheric scale height for the latitude inter- 
val 0° to -45° is shown in the lower panel. 

Similar to Figure 1 except for the latitude interval -30° to -40°. 

Similar to Figure 1 except for the latitude interval -40° to -50°. 

Amplitude and phase of depth averaged temperature employed as 
a zero order model. A phase of v/2 corresponds to a maximum 
at 18 hours local time. 

Amplitudes and phases of the zero and first order contributions 
to the meridional component of the diurnal wind field. The zero 
order contribution is a westward traveling wave while the first 
order contribution is an eastward traveling wave. 

Same as Figure 5 except for the zonal component of the diurnal 
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Figure 7 - Diurnal amplitude of meridional wind as a function of latitude and 
longitude. The amplitude shown results from a combination of the 
zero order and first order waves. The wind speeds shown are 
in m sec"^. 

Figure 8 - Same as Figure 7 except for zonal wind. 

Figure 9 - Vertical energy flux due to topographically excited vertically 
propagating tidal modes. 
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Figure 1. Comparison of atmospheric temperatures at the 2mb level with 
topographic relief. The temperature data points represent individual meas- 
urements obtained within the latitude interval -20“ to -30" for the local 
time intervals indicated. The solid lines represent the results of model 
calculations for -25" latitude and 6.5 and 17.5 hours local time. The mean 
surface, elevation expressed in units of atmospheric scale hieght for the 
latitude interval 0" to -45" is shown in the lower panel. 
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Figure 2. Similar to Figure 1 except for the latitude interval -30° to -40°. 
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Figure 4. Amplitude and phase of depth averaged temperature employed as a 
zero order model. A phase of n/2 corresponds to a maximum at 18 hours 
local time. 
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WIND AMPLITUDE (M SEC) PHASE (RADIANS) 


DIURNAL WIND 
MERIDIONAL COMPONENT 



LATITUDE 

Figure 5. Amplitudes and pliases of tlie zero and first order contril:)utions to 
tlie meridional component of the diurnal wind field. The zero order contri- 
bution is a westward traveling wave wWle the first order contribution is an 
eastward traveling wave. 
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LATITUDE 

Figure 6. Same as Figure 5 except for the zonal 
component of the diurnal wind. 
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